Nature of light correlations in ghost imaging 



Sammy Ragy and Gerardo Adesso 
School of Mathematical Sciences, The University of Nottingham, 
University Park, Nottingham NG7 2RD, United Kingdom 
(Dated: April 27, 2012) 

We investigate the nature of correlations in Gaussian light sources used for ghost 
imaging. We adopt methods from quantum information theory to distinguish genuinely 
quantum from classical correlations. Combining a microscopic analysis of speckle- 
speckle correlations with an effective coarse-grained description of the beams, we show 
that quantum correlations exist even in 'classical'-like thermal light sources, and appear 
relevant for the implementation of ghost imaging in the regime of low illumination. We 
further demonstrate that the total correlations in the thermal source beams effectively 
determine the quality of the imaging, as quantified by the signal-to-noise ratio. 




FIG. 1: Diagrammatic representation of a ghost imaging scheme. 



Ghost imaging [TJ [2] is an imaging modality configured as in figure [Tj Behind the object to be imaged, a 
bucket detector is placed. This detector yields no spatial information about the object but merely informs of 
the net photon count or intensity of the light impinging upon it. On its own it cannot provide an image of the 
object. In the second branch of the scheme, however, there is a spatially resolving detector (e.g. a CCD array); 
in this case, the light impinging upon this detector never passes through the object so it also cannot, on its 
own, yield an image of the object. In order to obtain an image, it is necessary to correlate the outputs of these 
detectors. Additionally, we need a source of correlated light to begin with, to provide the illumination; it is the 
exact nature of this source which has been the cause of controversy. 

The first demonstration of ghost imaging was carried out by Pittman et al. in 1995 [3J. In their paper, they 
used entangled beams produced by spontaneous parametric down conversion (SPDC) to perform the imaging. 
Even then, they noted briefly at the end of their paper that it may be possible to adapt their technique to 
'classical' light — a proposition which sparked debate When this was experimentally achieved by Bennink et 
al. in 2002 [5], part of the debate was resolved. Over the following years, even more theoretical and experimental 
papers ironed out the differences between ghost imaging with 'quantum' and 'classical' light 6 8 . Notably for 
our purposes, Gatti et al. [5] proposed thermal-light ghost imaging which was experimentally achieved shortly 
afterwards (TU1 HJ . This prompted a debate on the big question whether the physical origin of thermal-light 
ghost imaging can be explained entirely using classical intensity correlations and to this day disagreement still 
persists [121IT3]. 

Despite the large number of studies addressing practical questions such as signal-to-noise ratio (SNR), image 
contrast and acquisition time H1HTB] . to our knowledge, there has been a lack of study on the specific 
decomposition of classical and quantum correlations (in their broadest sense) in the light source used for ghost 
imaging. We believe that a study from such an angle will elucidate how quantumness versus classicality manifest 
in the scheme. 

One qualifier of 'classicality' for ghost imaging has been that it can be accounted for by a semiclassical 
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optics description — where light is treated classically, but the quantum nature of the detector is accounted for 
by consideration of the shot-noise [l4l [17] . Thermal-light ghost imaging can indeed afford such a description. 
However, an immediate issue arises: for two modes of Gaussian light between which any correlation exists at all, 
a nonvanishing part of the total correlations has necessarily a genuinely quantum nature [181 119| , as captured 
for instance by the quantum discord EI] ■ Since the correlated thermal light used for such experiments is 
indeed Gaussian, it warrants a study of the extent to which these quantum correlations feature in so-called 
'classical' ghost imaging. 

To put this on a more rigorous footing, in quantum optics, 'classical' light is taken to be light with a proper 
(non-negative and non-singular) Glauber P-representation, meaning it can accurately be represented by clas- 
sical optical theory [35]. For two-mode Gaussian states in standard form [23], classicality as qualified by the 
P-representation criterion coincides with separability, that is, absence of entanglement [24] . However, this clas- 
sification does not correspond to the most general informational definition for the classicality of correlations in 
bipartite states. According to a recently established paradigm, rooted in quantum foundations and information 
theory, a state is classically correlated when it has vanishing quantum discord ED] EI], i-e., when a projective 
measurement exists that once performed on a subsystem leaves the composite state undisturbed. In this respect, 
almost all unentangled states possess general quantum correlations as well [25] , Quantum correlations without 
and beyond entanglement are capturing the attention of a vast community of researchers, as those correlations 
have been linked to practical advantages in a number of noisy quantum communication and computation frame- 
works where entanglement is too fragile to be maintained |26[ 127] . In fact, on general grounds, a recent work 
has indicated that the sets of states corresponding to each definition of classicality (the optics one in terms of 
P-representation, and the informational one in terms of discord) are almost completely disjoint EE]. Since the 
ghost imaging scheme is entirely dependent on the correlations between each beam, one needs to apply special 
care when tagging any instance of that protocol as 'classical' based on one perspective only. Considering that 
ghost imaging is a scheme which entirely depends upon the correlations between each arm, it seems likely that 
an analysis from the informational perspective of correlations may shed some light on the problem. 

To better understand the debate, it is useful to divide thermal-light ghost imaging into two categories: lensed 
and lensless. At present, advocates for the quantum description of thermal- light ghost imaging propose that the 
ghost image formation in the lensless case cannot be explained semiclassically but can only be accounted for by 
non-local two-photon interference, a distinctly quantum phenomenon. It is the indistinguishability of different 
ways of triggering a joint detection which is crucial for this characterisation [12l [29] . More so, there have been 
claims that the non-local interference picture of ghost imaging is necessary to explain the results of some recent 
experiments. These include ghost imaging with speckle-free light [30 (although this idea of speckle-free light 
has met some resistance [13] ) and also an experiment on turbulence-immune ghost imaging [31] , but due to the 
general contentiousness of ghost imaging this too has come under scrutiny [131 132] . 

On the other hand, the lensed picture is much less controversial. Even in terms of the two-photon interference 
picture, it is considered as a "man-made" local correlation; the propagation characteristics induced by the lens 
(i.e. a Fourier transform), mean that each mode in the source plane is mapped to a unique position in the 
detection planes and a local, semiclassical picture can be used. A detailed description of this can be found in 
[33] . Nevertheless, this type of ghost imaging still warrants an examination from our perspective due to the 
aforementioned presence of quantum correlations in all correlated Gaussian states. 

In this paper, we quantitatively compare the correlations in the source Gaussian light for lensed ghost imaging 
to a specific figure of merit, the SNR. In the case of thermal- light ghost imaging, we find that there is a 
universal quasi-linear relation between the total correlations (the sum of quantum and classical) and the SNR. 
Furthermore, our analysis reveals two regimes: one — for low illumination — in which quantum effects dominate 
(in the sense that quantum discord takes up most of the total correlations) , and another for high illumination — 
in which it is the classical side of correlations which is much stronger. Increasing the illumination intensity we 
thus witness a quantum-to-classical transition in the description of ghost imaging, all the time using unentangled 
thermal light, which is classified as entirely 'classical' by optics standards. In the limit of very high illumination, 
thermal light and entangled sources yield the same performance for ghost imaging, and we find that they are 
characterised by the same, mostly classical correlations (although we re-emphasise that the quantum component 
of correlations is only exactly zero for completely uncorrelated light). We also examine the behaviour of 
correlations in near-field, lensless ghost imaging, removing any dependence on speckles from our calculations. 
We find that both lensed and lensless ghost imaging behave in a similar manner. 
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RESULTS 



Overview 



Before delving into our results, it is necessary to clarify which variables influence the light correlations and 
the SNR simultaneously for our lensed picture. To begin we have the speckle-count per pixel, M of the spatially 
resolving detector. When we refer to 'speckle', we mean the spatial intensity fluctuations projected upon the 
detectors by a spatially incoherent beam |34j . The second parameter of interest is the illumination / = M/i, 
where /i denotes the average photon count per mode — which for our lensed formulation corresponds to the 
photon count per speckle. The final parameter of interest is the pixel count of the image, R, which informs us 
over how many modes the bucket detector integrates. More details on this formulation of ghost imaging can be 
found in [35 , from which we have taken the SNR expression. Out of several possible choices [H3I3S], we choose 
to adopt the SNR corresponding to an imaging function which utilises the covariance of the intensities in each 
beam. The explicit expressions of the SNR for thermal-light and SPDC-entangled light sources are provided in 
the Methods section. 

In the lensed ghost imaging, we have pairwise correlated speckles. Each pair of correlated modes in the 
source-plane is mapped onto a unique position in each of the CCD and bucket planes. As such, any speckle is 
correlated to one — and only one — speckle on the opposite plane. This is a very typical setup for the observation 
of so-called 'classical' ghost imaging [TUl ITT1 l3"S"] . 





Superposed 

(d) 



FIG. 2: Conceptual diagrams indicating the areas we are considering on the CCD and bucket planes and a superposition 
of both. For part (a) we are considering exactly matched areas on either plane, for part (b) we are considering entirely 
disjoint areas. Part (c) has two areas of equal magnitude, half overlapping. The correlations between these are ~ those 
for an individual speckle within the overlapping area. Part (d) considers one area | the size of the other, completely 
enclosed by it on the superposed picture. This case notably, is most like the comparison of pixel and bucket detector. 
The effective modes also have \ the cross-correlations between them as for any individual speckles within the overlap 
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We shall look at the character of two different types of correlations: the individual (microscopic) speckle- 
speckle correlations for each pair of correlated speckles, and a 'coarse-grained' view which averages over the 
entire bucket detector for each pixel of the spatially resolving one. This latter view, which allows us to adopt an 
effective two-mode description of the problem, gives an impression of the correlations involved in shaping each 
pixel of the ghost image with a close regard to how the ghost imaging scheme is actually constructed. There 
are several note-worthy aspects of the coarse-grained formalism which we shall elucidate diagrammatically; a 
full description of the derivation can be found in the Methods section. 

We specialise here our qualitative picture to thermal-light lensed imaging. We define effective operators for an 
area on each plane, cla 1 and cla 2 - It turns out that for any given area, A, on either plane, the autocorrelations 
(a\a,A) are equal to those for any individual speckle within the area. If we look at any two matched areas 
on the CCD and bucket planes (figure [i|a)), then the cross-correlations, (a Ai aA 2 ), are again exactly equal to 
the cross-correlations of the individual paired speckles (all speckles are assumed to have identical statistics). 
Similarly if we look at any two unpaired areas, in figure [2jh) then we get no correlations at all, which is exactly 
as we would expect. In this sense the coarse-graining proves to be a sensible method of averaging. Better yet, if 
we look at partially overlapping regions and regions of disparate area, the cross-correlations depend on the ratio 
of the overlap to total area. It scales as ■ This accounts for the cases of figure 2^c) and (d). Panel (d) in 

/ a 

particular is of relevance for ghost imaging; in this case, we find that the cross-correlations scale as 



or equivalent^ ^3^ = ^. 

In the following, quantum correlations Q will be quantified by the Gaussian quantum discord [18H20] , classical 
correlations C by the complementary one-way correlation measure introduced in |21) . and total correlations T 
by the conventional quantum mutual information, which coincides with the sum of the two. 



Speckle-speckle correlations 

Looking at the speckle-speckle correlations, we get our first hints that even with non-entangled thermal light, 
the role of quantum correlations cannot be ignored. In figure |3j we clearly see that for low illumination / — or 
alternatively, high speckle-count per pixel M — quantum correlations can actually exceed classical ones. This 
stems from the fact that for individual pairs, the correlations only depend on the expected photon count \i = jj . 
When jit < 1 the speckle-speckle quantum correlations dominate over the classical ones, Q > C. Entanglement 
is never present in the considered light source, yet our study reveals a definite non-classical nature of such light, 
manifested in the correlations between individual pairs of speckle modes in the low illumination regime. 




(a) (b) 

FIG. 3: Speckle-speckle correlations in thermal-light lensed ghost imaging, plotted as a function of the illumination / 
and the speckle-count M. In (b) a detail at M=l is shown. Notice the intersection between classical and quantum 
correlations at M = I. Classical ones dominate for I > M, while quantum ones are relevant for I < M. 
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Coarse-grained correlations 

The results we find for the coarse-grained correlations are even more striking. In the effective two-mode 
description, we find it convenient, for the sake of the discussion below, to normalise the correlations multiplying 

them by a factor \J~^ (see Methods for a justification). In the following, we refer to the normalised total, 

quantum, and classical coarse-grained correlations by T, Q, and C, respectively. Our aim is to compare their 
behaviour with the SNR, which quantifies the quality of the imaging. In this way we can attempt to identify 
which aspects of the correlations in the source beams can be regarded as key resources for the protocol. 

To begin, we first plot the total normalised coarse-grained correlations T and the SNR on separate graphs in 
figure [4j Despite the very different (physical and mathematical) nature of the two quantities under scrutiny, it 
is immediately noticeable that they have a very similar form, a fact which is even more evident when we observe 
a parametric plot of the SNR versus T in figure Qc). By varying /, M, and R in their region of interest (in 
particular keeping R ^> 1, which means nontrivial imaging), we find that the SNR always exhibits a quasi- linear 




FIG. 4: Coarse-grained analysis for thermal-light ghost imaging (with pixel count R = 100). Panels [(a)-(b)]: Log-linear 
plots of (a) Signal-to-noise ratio (SNR) and (b) total correlations T for increasing illumination /. In panel (c) we display 
the tiny region filled by SNR versus T with varying M £ [1, 1000] and I 6 [0, 1000]; the quasi-linear region is essentially 
invariant upon variations of R as well. Panels [(d)-(e)] display the decomposition of coarse-grained correlations, for 
M = 1. Looking at (d) we can see that the total correlations are bounded, and at high illuminations they are mostly 
classical. Zooming in the low-illumination regime near the origin, though, we see in (e) that there is an interval where 
quantum correlations still exceed classical ones. All the correlations are normalised by \J R/2. 
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(a) (b) 

FIG. 5: (a) Coarse-grained analysis for ghost imaging using entangled light produced by SPDC (with pixel count 
R — 100): Plot of SNR versus total correlations for M = 1 and I £ [0, 1000]; compare with the corresponding thermal- 
light case, figure^c). Panel (b) depicts a comparison between total coarse-grained correlations T for the cases of thermal 
and entangled light sources (with R — 100) as a function of the illumination /; it is shown that they share a common limit 
in the regime of high illumination [see Eq. Q], although in this setup SPDC-entangled light is always more correlated 
than thermal-light. All the correlations are normalised by ^/R/2. 

dependence on the total correlations with slope ~ 1. The relation is not exactly linear, yet the discrepancy 
between SNR and normalised total correlations stays smaller than one percent in the relevant parameter regime. 
This indicates that for any of the parameters that affect the light correlations, the SNR is affected in exactly 
the same way. The small difference may reflect the fact that there are properties of the detector and object 
which affect the SNR, but do not influence the correlations at the source. Rigorously, recalling /i = I/M, and 
using the formulas provided in the Methods, we have 

SNR _^ ^1/(2^ + 1) _^ x (1) 
f m->°° y/R/2]n[R/(R-l)] 

Our analysis thus reveals that the joint contribution of quantum and classical correlations can actually be 
used as a predictor for the performance, as measured by the SNR, of ghost imaging with 'classical', thermal- light 
sources. It is now of interest to look at how the total correlations decompose specifically for our coarse-grained 
description. In figure |4|d), we see that for high values of illumination, the quantum correlations Q tend to 0. 
As such, total correlations are almost all classical (and, like the SNR, they are bounded from above). We do, 
however, still see a region near the origin [figure He)] where the coarse-grained quantum correlations exceed the 
classical ones, in agreement with the speckle-speckle analysis of figure [3] 



Extending the picture 

We have comprehensively unveiled the presence and role of genuinely quantum correlations in thermal- light, 
'classical'-like ghost imaging. We can also extend our study to the case of ghost imaging with entangled 
light sources produced by SPDC. In this case, it has already been pointed out in [TS] that the entanglement 
between the individual modes acts as an extra resource and can lead to an increase in the SNR compared 
to the thermal-light case. In our formalism, we find that for entangled light the SNR (see Methods) grows 
much faster than linearly as a function of the coarse-grained total correlations T, as shown in figure pla) . 
However, in the limit of very high illumination, the coarse-grained total correlations for both thermal-light and 
entangled cases converge to the same asymptotic maximum value, which is an intuitively expected result [14] . 
Similarly, the SNR for both cases converges to the same upper limit in the regime of high illumination, given 
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by SNR^oo = v / M/[6 + M(2i?+ 1)] [35]. Specifically, we find that Eq. Q holds exactly for entangled as 
well as thermal light. When we plot the (normalised, coarse-grained) total correlations for thermal-light versus 
SPDC-entangled sources on the same graph, as in figure |5jb) , it is evident though that in the case of entangled 
light the correlations are always higher for given finite values of the parameters /, M and R. 

So far we have made quantitative comparisons for lensed ghost imaging only, however it turns out that the 
correlations for near-field, lensless ghost imaging also scale in precisely the same way. Instead of considering 
speckle-count per pixel and pixel-count over the image, we simply calculate the correlations between pixel area 
and bucket detector area. The qualitative results of our previous description then still hold, and we have 
included detailed calculations in the Methods section. 

DISCUSSION 

We have analysed the nature of correlations in Gaussian light sources used for ghost imaging from a quantum 
informational perspective, combining a microscopic with an effective coarse-grained description. We have found 
that even so-called 'classical' thermal light contains nonzero genuinely quantum correlations — as measured by the 
quantum discord [HI [2D] — whose contribution to the performance of ghost imaging schemes has been assessed. 
Since the entire scheme is dependent upon the correlations between the arms, the quasi-linearity between the 
total correlations and SNR strongly suggests that the total correlations may be acting as the essential operative 
resource. 

The dominant strength of quantum correlations in low-illumination regimes has an immediate physical ex- 
planation. In the limit of low illumination, there are few photons per mode and this is a regime where the 
quantum behaviour of light becomes very apparent. In fact, in a recent review, the non-local description for 
Gaussian thermal light is formulated with the explicit assumption of low photon-flux [17] and in the original 
paper expounding the non-local picture [29] , the adopted model of thermal light implicitly assumes a photon- 
counting regime, which is an equivalent criterion to low photon-flux. As elucidated by the results of this paper, 
it is likely no mere coincidence that such an assumption need be made, but rather a consequence of the fact 
that the quantum component of correlations available for detection by our scheme vanishes in the limit of high 
illumination. 

This reveals an interesting feature associated to the coarse-grained formalism put forward here. It actually 
indicates how the quantum nature of the light source becomes quenched as we diverge from the photon-counting 
regime and enter the classical limit of intensity correlations. For these high illuminations, the quantum correla- 
tions available for detection by our scheme tend to zero, and the physical model of the scheme does not require 
a quantum description of the light to be accurate. The small-scale features of the speckle-speckle correlations 
cannot be detected by the way the light is averaged over the bucket detector. This is an unusually clear look 
into how a system can transition from a quantum regime into a classical one. 

We have briefly extended our analysis to the less controversial case of SPDC-entangled light sources. It 
has been observed in previous work |14j that for high brightness the results of the ghost image formed with 
'classical'-like thermal light are "excellent approximations for the quantum [entangled] case" and we find indeed 
that in this limit, the (coarse-grained) total correlations in the sources too approximate each other closely. 

Our calculations, as explicated in the Methods section, reveal furthermore that lensed and lensless ghost 
imaging behave in a similar manner to each other. In the recent debate, which is well characterised by a comment 
from Shih |12) with a rebuttal from Shapiro et al. |13| . we find our results support the latter's arguments that 
near-field lensless imaging can afford a classical interpretation. In high illuminations the quantum correlations 
are suppressed. Ultimately, this shows that the geometry of the system has a limited role in determining its 
classicality or quantumncss, whereas the illumination conditions are a clearly more significant factor. 

We finally mention that quantum discord has been very recently measured experimentally in Gaussian states 
of light |36 3l ;3S] . We believe our predictions on the role of quantum versus classical — and total — correlations in 
ghost imaging with thermal Gaussian light could be tested with current optical technology. 
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METHODS 



Correlations in Gaussian states 



A Gaussian state is any state which can be represented by a Gaussian Wigner distribution in phase space. 
Ranging from thermal to squeezed states, they are convenient not only for their experimental ubiquity and 
robustness, but also for their ease of use for calculations. 

Gaussian states are entirely characterised by their first and second statistical moments. For this paper, 
we deal with zero-mean states; as such, the relevant object for our computations is the two-mode covariance 
matrix given by cry = Tr [p(RiRj + RjRi)]. Here Rij are the elements of our quadrature vector given by 

R = (qi,Pi,q2,ih) T E3- 

From this, we can fully determine the information required for computing the correlations between two 
modes. The measure of quantum correlations which we choose is the quantum discord 20, 21 , which is defined 
for bipartite states as the difference between two classically identical definitions of mutual information. When 
the difference is non-zero, we then necessarily have quantum correlations in the state. 

Any two-mode covariance matrix can be written in the form 



fl2 



fa c 0\ 

a d 

c b 

\0 d bj 



Subsequently, we can write the quantum discord, the classical correlations and the mutual information, respec- 
tively, as 

Q{ei2) = f(b) - f(u+) - f{v-) + inf /(Vdrti), 

C(a 12 ) = f(a)-Mf(y/e), 

T(a 12 ) = f(a) + f(b)-f(u + )-f(v_), 

where f(x) = ln(^p) + \n(^^), v± are the symplectic eigenvalues of the covariance matrix [23] and e 
is the conditional covariance matrix of subsystem 1 after an optimal Gaussian measurement has been performed 
on subsystem 2. The full details on the interpretation of these quantities and on of how to perform these 
calculations can be found in |18j . 



Lensed Ghost Imaging and SNR 

For this section we have utilised the model of Brida et al. [35] and a detailed analysis can be found in 
their paper. In the lensed setup, modes with particular values of transverse vector k are projected onto unique 
transverse coordinates p in the imaging plane. We assume that the modes corresponding to each value of k have 
identical statistics. Then, for the thermal-light case we can model the speckled beam as a set of independent 
modes a^. Performing the beam splitter transformation then leads to new operators b\.i = (a, + a vac ,i)/V2 and 
t>2,i = ifli ~ a vaCi i) I 'v2- The statistics for each of these beams individually are identical. 

The SPDC-entangled case requires a very similar calculation. The difference is that our modes are generated 
by the well-known equations b\^ — f/di(c/,) -I- Va 2 {— qi) and 6 2 ,i = Ud 2 (qi) + Va\{— qi), where U 2 — V 2 = 1. 

It is now useful to consider the signal-to-noise ratio of ghost imaging. The SNR is defined as 

gjv^-j^ _ I (Sj n — S ou t) | 

v /var(S'. m - S out ) 

where our object is a binary amplitude mask; that is, we either have full transmission or full occlusion of the 
incoming light. Si n corresponds to being in the object profile (full transmission), and S ou t the opposite. Here 
S is whichever imaging function we choose. 

As mentioned in the main text, for our analysis we are interested in three parameters. The first is the number 
of speckles, M, per pixel of the spatially resolving detector. The second, is the illumination, / = Mfi, where p 
is our mean photon number per spatiotemporal mode, and the last is the number of pixels in our ghost image 
R. 
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Choosing to take our imaging function S to be the correlation function of intensity fluctuations (the covariance) 
we have S(p) = (Ni — (Ni))(N 2 (p) — (N 2 (p))). The spatial dependence of N 2 refers to the position in the 
spatially-resolving arm. Furthermore, Ni denotes the net photon count on the bucket detector. To clarify, we 
have Nj = Ylf=i ^l,j> where Nij, refers to our photon count on the j th spatial resolution cell. In turn we can 
decompose each of our Ni, N 2 into the sum of individual modes (speckles) such that Ni — Ylj=i n i,j where 
i = 1,2. This quantity tells us the number of speckles per pixel. 

We then have [35] 

SNRthcrmal = — , A* ' 

\J (m 2 (2MR + M + 6) + 4ft(MR + 1) + 2MR + 1 

SNR rn t an gipH = — , , (3) 

s ^p 2 (2MR + M + 6) + J42MR + M + 6) + 1 V ; 

for thermal-light and SPDC-entangled light sources, respectively. 



Coarse-grained description for lensed imaging 

As outlined earlier, the SNR depends on an imaging function of the form S(x) = (Ni — (Ni))(N 2 (x) — (N 2 (x))), 
which takes into account multiple modes. Whilst it is easy to generate multimode covariance matrices, the 
problem of calculating the discord from these matrices is quite intractable. As such, we proceed to define an 
effective two-mode matrix. The key to this is to produce averaged operators which, under the expectation value, 
behave like usual single-mode operators. 

To wit, starting with a pixel which collects M modes and noting that N t — J2^Li n i,j — (f° r 
« = 1, 2), we wish to find mode operators c, for each pixel, such that {c\cj) — (2j=i ln order to do this, 

we need to explicitly consider the spatial dependence of the modes at the detection planes. Doing so, we can 
say Cj = YljLi i>i,jfi{p — Pj) where pj corresponds to the position of the j th mode on the orthogonal plane. Note 
that we do not consider the effects of the transmission mask on light propagation in this particular calculation. 

These operators then behave as we desire, but for one very crucial part: the commutation relations. This 
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Ci, then we 



is easily remedied though. If we define our effective operators in the following manner: di 

find ([d, d']) = 1. Performing similar steps for pixels collecting R modes, and then transforming our covariance 
matrices from the mode operator to the quadrature basis, enables us to obtain the effective 'coarse-grained' 
covariance matrices given by 



in 




/1 + 2/i 
l + 2p 



\ 

2^ 





and 



l + 2p 







Vb. 




l + 2p 

1 + 2p 













1 + 2/i } 



(4) 



for thermal-light and SPDC-entangled light cases, respectively. 

Notice that the off-diagonal blocks of the covariance matrices, which encode the intermodal correlations, 
appear scaled by a factor oc 1/ y/~R as a consequence of the averaging applied to preserve the commutation 
relations. This makes e.g. the covariance (dJ A &a 2 ) sensitive to changes in one of the areas, which should not 
be, as an increase in, say, A 2 for a fixed A\ < A 2 implies no loss of any physical correlation. Therefore, to make 
sure we measure genuine phase-insensitive correlations between the beams in the effective two-mode description, 
we can a posteriori renormalise classical, quantum, and total correlations, evaluated on the matrices of Eq. Q, 
multiplying them by a factor proportional to V^R. The specific choice \J R/2 is dictated by mere convenience, as 
it makes the slope in the SNR versus total correlations relation converge exactly to 1 [see Fig. |4|c) and Eq. Q]. 
Needless to say, the quasi-linear interdependence existing between the two quantities for thermal light, as well as 
all the results derived in the main text and illustrated in figures [4] and [5] are obviously not qualitatively affected 
(albeit for quantitative rescaling) by any specific normalisation procedure implemented on the correlations. 
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Coarse-grained description for lensless ghost imaging 



The above formalism for lensed ghost imaging considers a simple one-to-one mapping from the operators in 
the momentum-basis at the source plane, to the position dependent ones at the detection planes. However, we 
can generalise our calculations to other sorts of propagation. For the purposes of lensless ghost imaging, it is 
conventional to consider free-space propagation in the paraxial approximation. In this case, it is possible to 
characterise the field at transverse point pj in plane Zj by E(pj, zf, k) cx J2 k 9{pji z ji k)ci(k) for j — 1, 2, where 
cl(k) are our source-plane operators and g(fij, Zf, k) is the Green's function describing the propagation of the 
field to the point with transverse coordinates pj on the detection plane |33[ 139] . The form of g(pj, Zj] k) in the 
quasimonochromatic case is [39] : 



— ihr,p ik a z i f - As i |2 

9{ Pj ,z r ,K) = J dp s e^- p ' 1 e-***. 

We can then calculate the the auto-correlations and cross-correlations and from this, our effective covariance 
matrix. To begin, we note that in the quasimonochromatic, paraxial approximation, [Ei(pi, Zj), Ej(pj, Zj)] = 
6(pi — Pj)Sij, where p are the coordinates in the transverse plane [H] HO]. As such, these behave under the 
expectation value as though they were the usual single mode operators a. Henceforth, we will assume that we 
are always considering correlations on planes at equal distances from the source as is conventional for ghost 
imaging setups. For individual modes or correlated pairs: 



{E\p h z)E{ Pj ,z))K 



^ 9*(Pi,z; K)g{ Pj ,z; n')a) (k)ci(k')) 

kk' 



k e 



ik()Z 



kg 
2nz 

fcp 
2?rz 



2ttz 

2 



dp s e l 2 



\Pi-Ps\ pin-Pa 



dp' s e l ^~ p '^ e- iK -»'°)(at(K)a(K)) 



dp s dp' s e 1 ^ 



> p iH\Pi~P' s \ 2 



~\Pi-Ps\ gift'Cpa 



-P's)\ 



(a t (K)a(K))( / dp s e 



(J^j 2 (at( K )a( K ))e^(P?~P?) J dpj^"^-^ 



2nz 

(a* (K)a(n))S(p 



For the calculation, we have assumed a large disk-like source and a large number of modes. These are standard 
assumptions made in the derivation of the non-local biphoton model of ghost imaging 29, 33]. Our calculation 
indicates that (pi, z)E(pj, z)) cx (a) (k)ci{h)) when pi — pj, and is otherwise. The correlation at any point 
on a single plane or paired points on CCD and bucket planes is proportional to the expectation value for a 
given mode in the momentum-basis. This behaviour reflects that of lens-based ghost imaging. In fact, it shows 
that we have a proportionality between photon counts in the source-plane and detection planes, which is not 
an entirely surprising result. 

We can also go further and define effective operators in a similar manner to how we did for lens-based imaging 
by writing E p = I — J A E(p, z)dp where we are integrating over a pixel with area A p . We can test that the 



correct commutation relation holds for this operator under the expectation value using similar assumptions as 
above. 

It is then also easy to show that (E^E p ) cx (a' (k)ci(k)) . Again, the behaviour reflects that of lensed ghost 
imaging. This leaves an integration over the bucket detector to fully characterise our scheme. The only difference 
to the above calculation is that we are integrating over a larger area, A\,. It is evident that the auto-correlations 



will be the same. For the cross-correlations we find {E h E- 



Essentially, we scale the cross- 



correlations by the inverse square root of the ratio of the bucket detector area to the pixel area. This yields an 
identical scaling with changes in our three parameters /, R, M as for the lens-based imaging. 
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